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We present a detailed derivation and testing of our approach to rescale the dynamics of mesoscale 
simulations of coarse-grained polymer melts (I. Y. Lyubimov et al. J. Chem. Phys. 132, 11876, 
2010). Starting from the first-principle Liouville equation and applying the Mori-Zwanzig projec- 
tion operator technique, we derive the Generalized Langevin Equations (GLE) for the coarse-grained 
representations of the liquid. The chosen slow variables in the projection operators define the length 
scale of coarse graining. Each polymer is represented at two levels of coarse-graining: monomeric 
as a bead-and-spring model and molecular as a soft-colloid. In the long-time regime where the 
center-of-mass follows Brownian motion and the internal dynamics is completely relaxed, the two 
descriptions must be equivalent. By enforcing this formal relation we derive from the GLEs the 
analytical rescaling factors to be applied to dynamical data in the coarse-grained representation to 
recover the monomeric description. Change in entropy and change in friction are the two correc- 
tions to be accounted for to compensate the effects of coarse-graining on the polymer dynamics. 
The solution of the memory functions in the coarse-grained representations provides the dynamical 
rescaling of the friction coefficient. The calculation of the internal degrees of freedom provides the 
correction of the change in entropy due to coarse-graining. The resulting rescaling formalism is a 
function of the coarse-grained model and thermodynamic parameters of the system simulated. The 
rescaled dynamics obtained from mesoscale simulations of polyethylene, represented as soft colloidal 
particles, by applying our rescaling approach shows a good agreement with data of translational 
diffusion measured experimentally and from simulations. The proposed method is used to predict 
self-diffusion coefficients of new polyethylene samples. 



I. INTRODUCTION 

The past few years have witnessed a growing interest 
in the design and application of coarse-graining meth- 
ods to simulate complex fluids. [T] This effort has been 
motivated by the need for improving computational ef- 
ficiency with the purpose of investigating complex sys- 
tems on the numerous lengthscales on which their prop- 
erties develop. I2HS] Computer simulations have the ca- 
pability of providing detailed microscopic information on 
the static and dynamics of the systems under study, [7 
but they are limited in the range of timescales and in the 
number of molecules that can be simulated because the 
precision of the calculations degrades with the number of 
computer iterations with a behavior that depends on the 
Lyupanov exponent of the system. Once the number of 
particles is set, the window of achievable timescales that 
can be investigated becomes defined. [HI [H] Because the 
maximum number of iterations decreases with increasing 
number of simulated particles, it is particularly difficult 
to simulate systems where characteristic lengthscales are 
diverging, such as a system approaching a second order 
phase transition. [TOl [11] 

Recent improvements of computational machines has 
lead to a considerable extension of the maximum time- 
and length scales that can be reached by simulations 
where the system is described at the atomistic level. 
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However, for many complex systems, including liquids 
of high-molecular weight macromolecules, the computa- 
tional power is still inadequate to describe, at the atom- 
istic level, the long-time dynamics. For example, the 
most recent and advanced simulations of long chains that 
have an extended number of entanglements adopt a sim- 
plified model, which treats the structure of the polymer 
as a collection of beads and springs interacting through 
a FENE potential. This model allows for the simulations 
of a large number of polymers, which is important for the 
proper calculation of viscoelastic properties, and reaches 
full relaxation for all but the longest chains simulated. [1^ 

Progress has been made when the focus is on qualita- 
tive behavior and scaling exponents. [TUl [T5] For example, 
if the complex intra- and inter-molecular non-bonded in- 
teractions are simplified into an identical potential, the 
computational efficiency improves dramatically as the 
code does not need to identify and treat uniquely differ- 
ent pairs of interacting sites. This strategy, however, has 
the disadvantage that the thermodynamics of the system 
is not properly described because the interactions are too 
drastically simplified. 

The need for methods that are fully predictive of the 
physical properties of a system on the basis of the specific 
chemical structure of the sample and its thermodynamic 
conditions has stimulated new interest in developing fast 
quantitative simulations. Such predictive approaches are 
useful, for example, to evaluate a priori the structure 
and dynamics of newly synthesized polymeric materials, 
in relation to their technological applications. Following 
this perspective, several procedures have been proposed 
to speed up atomistic simulations, while conserving their 
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power of predicting quantitative properties. |141 A few 
simulations of long entangled chains have been performed 
using united atoms fUAl . [Til [TMT5] . For UA the effec- 
tive unit is very close to the atom is size, i.e. CH^ with 
X = 1,2,3, which allows for some gain in the computa- 
tional time. 

A useful strategy to improve the outcome of simu- 
lations on the long timescale and large lengthscale is 
the use of coarse-graining procedures. [21 [19] A coarse- 
graining procedure averages out irrelevant degrees of free- 
dom, which occur on lengthscales smaller than a desig- 
nated cutoff length, and this allows for the extension to- 
wards large scales of the simulations. Another way to 
put it is that, because the interaction potentials become 
softer, the maximum time and lengthscale increase as 
the basic timestep of the mesoscale (MS) simulations be- 
comes larger. The characteristic lengthscale of coarse- 
graining has to be defined on the basis of the properties 
that need to be investigated. In this paper we discuss a 
first-principles way of selecting meaningful lengthscales 
for the structural and dynamical coarse-graining. 

Several considerations need to be made to properly de- 
velop a coarse-graining procedure. As the coarse-grained 
liquid is represented as a function of new coordinates, 
an effective potential needs to be derived to be used as 
an input to the MS simulation. Care has to be taken to 
make the potential reproduce the structure of the system, 
namely pair distribution functions, and to be thermody- 
namically consistent. A common procedure to optimize 
the coarse-grained description is to use self-consistent nu- 
merical methods that are optimized to reproduce atom- 
istic descriptions through iterative procedures. Usually 
the target is the optimization of specific quantities, such 
as the pair distribution function, [2U] the forces generated 
by the soft potential, f6^ or directly the thermodynamic 
properties. 2U 

We recently proposed an approach that starts from 
the Ornstein-Zernike (OZ) equation where the atomistic 
sites are defined as real sites, and the coarse-grained 
sites are defined as auxiliary sites. [22H25] Because our 
procedure is analytical, and no optimization of param- 
eters is needed in our approach as the potential is ex- 
plicitly dependent on the thermodynamic and molecu- 
lar parameters, it opens up the possibility of deriving 
a formal solution to key problems. For example, it is 
straightforward to show that the structural properties 
are consistent between the two levels of description, i.e. 
atomistic and coarse-grained. [53j Moreover, the thermo- 
dynamic properties of the coarse-grained polymer liq- 
uid (e.g. isothermal compressibility. [22j pressure in the 
virial and in the compressibility routes, total and cohe- 
sive energy) are shown to be formally consistent in the 
two levels of coarse-graining. [29] Local structure is eas- 
ily included a posteriori through a multiscale modeling 
procedure. [551 1301 Finally, it is possible to derive an an- 
alytical rescaling factor for the dynamics, which is the 
main focus of this paper. 

While the structure is well described by simulations 



of the coarse-grained system on the scale larger than 
the scale of coarse-graining, the dynamics in MS sim- 
ulations is unrealistically fast. Because local degrees of 
freedom are averaged out, the coarse-grained molecules 
move rapidly over a simplified free energy landscape. As 
the system explores efficiently this "reduced" configura- 
tional landscape, the measured dynamics is artificially 
sped up by the smoothness of the potential. This is useful 
when coarse-grained representations are used to rapidly 
reach an equilibrated state of the system before start- 
ing the atomistic molecular dynamics (MD) simulation. 
However, to directly collect information on the dynam- 
ics of systems from MS MD simulations, it is necessary 
to develop formalisms that rescale the unrealistically fast 
dynamics into the slower dynamics at atomistic resolu- 
tion. In this paper we discuss in details an analytical 
procedure we recently proposed to rescale the mesoscale 
dynamics. The procedure is able to predict center-of- 
mass dynamics in quantitative agreement with experi- 
ments and atomistic simulations. 31J 

The common strategy to rescale the dynamics is to 
build a "calibration curve". The latter is obtained 
through the numerical fitting of dynamical quantities and 
optimization of the related parameters until the agree- 
ment of dynamical properties calculated in an all-atom 
and in a MS simulations is obtained. [32 133] However, 
the numerical calculation of optimized calibration curves 
for the dynamics is quite difficult to achieve for macro- 
molecular systems, as the dynamics is mode dependent: 
there are in principle N internal modes in any molecule 
formed by N units and the degree of polymerization of a 
long chain can be of the order of one million monomers. 
Moreover, numerically optimized parametric quantities 
are in general not transferable between systems in dif- 
ferent thermodynamic conditions or with different chem- 
ical structure or increasing degree of polymerization. To 
overcome this problem, it is common to select as coarse- 
grained units ones that are very close in size to the atom- 
istic units, so that the needed corrections to reach con- 
sistency in dynamic properties are minimal. In this case, 
corrections to the measured dynamics can be evaluated 
through a perturbative formalism which should rapidly 
converge to the desired value. The resulting computa- 
tional gain is, however, still limited. Recently a numer- 
ical Ornstein-Zernike-based approach, with atomistic- 
level coarse-graining and the Percus -Yevick closure ap- 
proximation, has been proposed, which shows different 
rescaling factors depending on the time correlation func- 
tion under study. [31] Another coarse-grained approach 
for polyethylene melts describes a polymer chain as a 
collection of soft blobs connected by elastic bands, which 
enforce chain-chain uncrossability. Simulations follow an 
effective Langevin equation, whose parameters, i.e. ef- 
fective potential, frictions and random forces, are ob- 
tained from an atomistic MD simulation. The optimized 
equation of motion (eom) reproduces well experimental 
data. [35] 

Our approach is different from others in several ways. 
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First of all it is analytical rather than numerical, pro- 
viding the formal rescaling factor by solving the eoms 
in the two levels of representation. In this way, there is 
no need of performing an atomistic simulation to input 
numerical quantities in our formalism. The dynamics 
measured in MS simulations of coarse-grained systems 
is directly rescaled into its atomistic counterpart using 
approximate closed-form expressions of friction and en- 
tropy. The two levels of coarse-graining, which allow for 
a straightforward analytical solution, are here two simple 
isotropic models: a soft sphere description for the molec- 
ular coarse-graining, and a bead-and-spring description 
for the monomer level coarse-graining. More sophisti- 
cated coarse-grained models can be developed for in- 
termediate lengthscalesP^ \Tl\: however, the formalism 
can become more involved. At the atomistic level the 
polymer is described as a "bead-and-spring" type of ap- 
proach where the chain is a collection of friction points 
connected by harmonic springs. This is an implemen- 
tation of the most popular model to treat unentangled 
polymer melt dynamics, i.e. the Rouse model, and maps 
well into the dynamics of polymers described not only by 
UA simulations, but also by atomistic simulations and 
experiments, as it contains both local chemical structure, 
semiflexibility, and finite size effects. [37H39] It is a very 
accurate and molecular specific model, which has been 
shown to describe well, for example, the dynamics of the 
protein CheY, by testing its predictions of NMR relax- 
ation against experiments. [37l l40] 

In our coarse-grained model a polymer chain is repre- 
sented as a soft-colloidal particle. p2 l I42ti44] Because the 
lengthscale of the coarse-graining is of the order of the 
molecular radius of gyration, i.e. the size of the molecule, 
the direct predictions of the rescaling procedure are suit- 
able for properties on lengthscales larger than Rg and 
on timescales longer than the longest time of intramolec- 
ular relaxation, i.e. the longest correlation time in the 
Rouse theory. Internal dynamics cannot be obtained di- 
rectly from the coarse-grained simulation; however, the 
rescaled diffusion coefRcient leads to the monomer fric- 
tion coefficient, which can be used as an input to well- 
tested theories of polymer dynamics, and indirectly re- 
covering the dynamics in the complete spectrum of poly- 
mer relaxation. An example of this kind of calculation 



is presented in this paper in Section VI C The extended 



level of coarse-graining provides a good test of our pro- 
cedure, as large deviations could result from the rescal- 
ing if the method were not correct. Furthermore, our 
procedure can be useful in the study of long-time relax- 
ation, given that large length scales and long timescales 
are most difficult to simulate for polymeric systems. 

Although the outline of our rescaling theory has 
been published recently in a short paper, [31 this paper 
presents a detailed derivation and discussion of our ap- 
proach, which includes the prediction of the dynamics 
for new samples. After introducing our coarse-grained 
model input to the mesoscale simulations, we formally 
derive the rescaling approach for the dynamics, starting 



from the Liouville equation and using projection oper- 
ators. Friction coefficients in the two descriptions are 
derived from the solution of the memory functions, while 
the rescaling of the simulation time is obtained from the 
entropic contribution, which accounts for the intramolec- 
ular degrees of freedom neglected in the soft colloid repre- 
sentation. Theoretical predictions compare well against 
UA MD simulations, [14. .45. -47] and experiments, |39, 4S1- 
[5T] . We also calculate the diffusion coefficient for new PE 
samples in thermodynamics conditions for which UA-MD 
data are not available. The purpose of these calculations 
is to show that our method is not a simple rescaling of 
the mesoscale data through a shift of the diffusion coef- 
ficient to bring dynamical results to coincide with atom- 
istic simulations, as is conventionally done. Instead our 
approach is fully predictive and can be used to calcu- 
lated the diffusion coefficient, and the monomer friction, 
for new samples. 

The paper is structured as follows: after introducing 
our coarse-grained model in Section [ll] and the projec- 
tion operator technique to derive the equations of mo- 
tion in the two levels of coarse graining in Section [TTlj we 
formally derive the rescaling approach for the dynam- 
ics from the solution of the memory functions in Section 



IV We then present the MS simulations (Section |V|), as 
well as the results obtained from the same, and apply 
the rescaling procedure to the data from MS simulations 



(Section VI I. Predictions of dynamical quantities and di- 



rect comparison for several samples, both from atomistic 
simulations and from experiments, provide a stringent 
test of the approach and show good quantitative agree- 
ment. A brief discussion in Section IVIII concludes the 
paper. 



II. COARSE-GRAINING OF POLYMERIC 
LIQUIDS: STRUCTURAL PROPERTIES 

In this section we briefly review the theoretical back- 
ground of the pair distribution functions that are in- 
put to our rescaling equation. The structure of a poly- 
meric liquid, at lengthscales equal or larger than the 
monomer size, is fully specified by the momomer total 
distribution function, h(r), which for polymer melts de- 
pends on two characteristic lengthscales, namely the den- 
sity fluctuation lenthscale, which is the atomic length, 
and De Gennes' correlation hole lengthscale. [52^ The lat- 
ter is o f the order of the molecular radius-of-gyration, 
Rg = y/N/6l, which is the overall dimension of the poly- 
mer, where N is the degree of polymerization and I is 
the statistical bond length. We select I and Rg because 
these are the two lengthscales that define the structural 
properties of the polymeric liquid. 

At the monomer level traditional dynamical ap- 
proaches, such as the Rouse model and semiflexible mod- 
els, adopt a bead-and-spring representation where each 
monomer can be modeled as a friction point connected 
by springs (see Fig. [T]) . A similar model, where the poly- 
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FIG. 1: Illustration of monomer and overall coarse-graining 
of a homopolymer linear chain. 



meric chain is described as a collection of "sites" centered 
at the center of the monomeric unit, is also in conven- 
tional theories of polymer liquids. [551 IM] Although "site" 
is the word most used in the liquid state community and 
"monomer" or "bead" is the common wording in the lit- 
erature on polymer dynamics, in this paper they identify 
the same CHx unit and henceforth they will be used in- 
terchangeably. It is important to notice that all the CH2 
units are assumed to be equivalent and independent of 
the position along the chain. 

The coarse-graining of a polymer at the Rg lengthscale 
represents each molecule as an interacting soft colloidal 
particle with symmetric or asymmetric shape. [1^44j In 
our model. [22H25] the macromolecular liquid is repre- 
sented as a liquid of symmetric soft colloidal particles 
interacting through a pair potential. This potential has 
a range of the order of few Rg, and each soft-colloidal 
particle is centered at the center-of-mass of a polymer 
(see Fig.[l]). 

The coarse-graining procedure that translates the 
monomer description into the solf-colloidal representa- 
tion is performed starting from an Ornstein-Zernike equa- 
tion where monomers are assumed to be real sites, while 
the center-of-mass (cm) are auxiliary sites. [JT] The cm- 
cm total intermolecular correlation function is expressed 
as a function of the polymer parameters as [22 [25] 




where erfc(a;) is the complementary error function. Here 
C = ^s/(27rp:J = 3/{7TpP) with = being 
the reduced molecular number density, pch = p/N is the 
molecular density, p the site number density, and I is the 
statistical segment length. The length scale of density 
fluctuations, ^p, is defined as = +?p~^i and = 
Rg/\/2 is the length scale of the correlation hole. [52] 



In the limit of long chains, N ^ 00, Eq.Q reduces to 
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For polymers with N > 30, a plot of h{r) shows that the 
two equations, Eqs.Q and ([2]), are indistinguishable. \22\ 

1131 

The structure of the liquid on the lengthscale of the 
polymer radius-of-gyration and larger, as represented by 
h'^'^{r), is in quantitative agreement with the output of 
both the atomistic UA MD and the MS MD simulation 
of the coarse-grained liquid. The theory recovers iden- 
tical analytical expressions of the compressibility in the 
atomistic and the coarse-grained representations, indicat- 
ing thermodynamic consistency between the two levels of 
description. [221 [25] 

Eqs.Q and ^ are de facto coarse-graining equations, 
which translate the atomistic description of a polymer 
liquid, onto its representation as a liquid of interact- 
ing soft colloidal particles of size Rg. The advantage of 
our coarse-graining approach is that it is analytical and 
general as it applies to systems with different thermody- 
namic conditions, different degree of polymerization and 
different bond length. 



III. 



DYNAMICAL COARSE-GRAINING: FROM 
THE LIOUVILLE TO THE LANGEVIN 
EQUATIONS 



While the structure of the polymeric liquid, as repre- 
sented by the total correlation function, is identical in 
the atomistic and coarse-grained descriptions, [22, 23] the 
dynamics of the coarse-grained system, as measured in 
the MS MD simulations of the soft-colloidal particles, 
is unrealistically accelerated. In Fig. (Ill I we show, for 
a polyethylene chain with N = 44, the mean-square- 
displacement of the center-of-mass obtained in MS MD 
simulations of the polymer liquid represented as soft- 
colloidal particles and the mean-square-displacement di- 
rectly measured in UA MD simulations . The dynam- 
ics in the coarse-grained representation is several orders 
of magnitude faster than the atomistic description. Be- 
cause the level of coarse-graining of the model presented 
here is extended, this effect is more evident than in other 
models; however, accelerated dynamics is present in any 
simulation of coarse-grained systems. 

It has been argued that there are two main effects of 
coarse-graining that accelerate the dynamics: namely, 
the change in entropy and the change in the friction 
coefhcient. Ottinger has presented an approach for 
systems far from equilibrium that accounts for those 
effects. [55] We propose here a procedure, based on first- 
principles theory to properly account for both contribu- 
tions through the introduction of the necessary correc- 
tions for systems where the fluctuation dissipation theo- 
rem applies, e.g. close to equilibrium. 
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ration of timescales is observed, it is necessary to in- 
clude corrections to the projected dynamics, which ap- 
pear as contributions to the friction coefficient, expressed 
as memory functions. In the system investigated here, 
polymer melt dynamics, no clear separation of timescales 
occurs between the dynamics of the "tagged" chain and 
the dynamics of the surrounding molecules. |59j For this 
reason, the Generalized Langevin Equation generated 
from this procedure needs to account for the correction 
terms to the projected dynamics, which are represented 
by the memory function contributions . [55] 

For a liquid of n macromolecules containing N 
monomers, the first-principle Liouville equation is sim- 
ply written as 



df (R,P,^) 

dt 



= ^L/(R,P,^) 



(3) 



with 



FIG. 2: Cm mean-square displacement, for a polyethylene 
melt with N = 44, from MS MD simulations (triangles) and 
UA MD simulations (squares). Dashed lines show the diffu- 
sive limits of the two samples. 



To coarse-grain the dynamics of the polymeric liquid 
on the lengthscale of the radius of gyration, we adopt a 
Mori-Zwanzig projection operator technique, where the 
selected slow variables are the position and momentum 
coordinates of the polymer center-of-mass. This descrip- 
tion should represent well center-of-mass diffusion. [56J- 
[58] 

The atomistic level representation is obtained follow- 
ing the same Mori-Zwanzing procedure, but choosing 
as the slow relevant variables the ensemble of position 
and momentum coordinates of the center of mass of 
the monomeric unit, which for a polyolefin is the CHx 
unit, with x = 1,2, or 3. This model is consistent with 
the representation of the polyethylene chain in UA MD 
simulations. [5^ [SO] and it has been shown to describe at 
a high level of accuracy the dynamics of polyolefins at 
the monomer lenghtscale. [STj [62l |63] 

In the long-time regime the two descriptions, soft- 
colloid and monomeric /U A, should be identical as they 
both recover the diffusive dynamics of the center-of- 
mass. [64 In fact, they are not, as the soft-colloidal de- 
scription is heavily coarse-grained and its dynamics is 
accelerated. The analytical rescaling factor is derived di- 
rectly from the comparison between the soft-colloid and 
the monomer dynamical equations. 

As this coarse-graining and rescaling procedure is gen- 
eral, it can be adopted to formalize the dynamics of the 
molecular liquid at the desired level of coarse-graining. 
However, the projection operator technique rests on a 
separation of timescales between the slow relevant vari- 
ables onto which the dynamics is projected, and the fast 
irrelevant variables that are averaged out. If no sepa- 



/(R,P,t) 



n 



N 



\{5{vi{t)-WMl>i{t)~Vi) 



,(4) 



the instantaneous distribution in reduced phase space, 
and RJj and are the phase-space variables associated 
with the Cartesian position and momentum coordinates 
of the bead a belonging to molecule i, namely v\{t) and 
pjj(i). The formal solution of Eq.(|3| is 



/(R,P,t)=e-*^7(R,P) , 

with the shorthand notation / (R, P) = / (R, P, 0) 
The Liouville operator is defined as 



(5) 



N r 



EE 



d 
dpi 



Pi 



d 
dri 



(6) 



where the total energy Uj in the Hamiltonian, con- 
tains both intramolecular, Uj, and intermolecular, Wij, 
pairwise decomposable potential contributions. The in- 
termolecular potential contains both interactions be- 
tween the n tagged chains, Wj^., and between the tagged 
chains and the surrounding ones, Wjk, so that the usual 
condition applies that Lof{Il,P) = 0. The statistical 
average of the phase space density is defined as 



/ (R, P)) = JdrJ dpf (R, P) (r, p) 



(7) 



with the equilibrium distribution of particle positions and 
coordinates 



dr / dpe 



-flH 



(8) 



where /3 — (ksT)'^, the Boltzman's constant, and 
T the absolute temperature. Following Mori-Zwanzig, 
we define the projection operator, P, for the coarse- 
grained model we adopt, namely the monomer and the 
soft-colloidal. 
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Monomer level representation of the polymer 
chain 



In our atomic-level description each macromolecule is 
represented as a collection of connected beads, or friction 
points. In the field variables for one molecule {n — 1), 



g{-R,-p,t) 



N 



(9) 



the projection operator is defined as 

Ph (R,P,<) = J dR' J dP' J dR" J dP" 

{h{R,P,t)g{R',P')) 
x{g{R',P')g{R",P")y'g{R",P") , 



(10) 



where P = (P)^ and Pg (R, P) = g (R, P). Here we use 
for the field variable the symbol g{R,P,t) to indicate 
that the slow variables in the projection operator can 
be different than the ones in the general formalism of the 
preceding section. By applying the projection operator to 
both the left and the right sides of the Liouville equation, 
one recovers a generalized Langevin eauation. |57H60j 

Briefly, the generalized Langevin equation in the phase 
space is then transformed into its analog equation in 
space coordinates, yielding 



' dt^ 



= f3- 



d 



dVa{t) 



ln-0(r) 



t N n 



b=l 



F?(i~r)) 



(11) 



where ip{r{t)) is the intramolecular distribution function. 
The inertial contribution in Eq.( 11 ) can be discarded, as 



the liquid has a low Reynolds number and the dynamics 
is overdamped. The Generalized Langevin Equation is 
simply written as 



Cr 



dVait) 



1 d 



ln'0(r) 



(12) 



dt 13 dVa{t) 

with the averaged friction coefficient, in the Markov limit. 



N 

E 

a,b=l 



dr(F,(t).F^(t-r)) . (13) 



This equation describes how the monomer friction coef- 
ficient is generated from the space and time correlation 
of the random forces that act on two different segments 
of the "tagged" polymer chain, a and h. The extent of 
the correlation depends on the propagation of the forces 
through the macromolecule, its structure and local flex- 
ibility. The forces are generated by the monomers of 
the surrounding molecules randomly colliding with the 



monomers of the tagged chain: the collision strength de- 
pends on the structure of the liquid and on the interpar- 
ticle potential. A more explicit definition of the friction 
coefficient is given in the following sections. 



B. Solution of the Generalized Langevin Equation 
in the monomer representation 

The intramolecular distribution function is approxi- 
mated in our description by a Gaussian distribution 



V'(r) = [(27r)^det(A-i)] e 



rr" Ar 



(14) 



which holds for polymer chains longer than about 30 
monomers. [37] This leads to a Generalized Langevin 
Equation where the intramolecular contribution is linear 
in the monomer coordinates 



Cr 



dVajt) 
' dt 



Ot. rp N 

-^^A,,,r,(t)+F?(0 , (15) 



b=l 



and is simply solved through transformation into normal 
modes of motion. [37l [64] The matrix A is defined, for 
a semiflexible polymer represented as a Freely Rotating 
Chain (FRG), as the product of two matrices, M and U, 



M 



(16) 



with the connectivity matrix, with dimensions N x N, 
defined as 



M 







... iv- 


'\ 


-1 


1 


... 





-1 1 


... 




V o' 





-1 1 


) 



(17) 



and the U matrix defined as a function of the stiffness 
parameter g as 



U, 




(18) 



Here, g — —(cos 6*) and 9 is the angle between 
two consecutive bonds in the FRC representation of a 
homopolymer. |64j The stiffness parameter, g, is specific 
of the chemical structure and thermodynamic conditions 
of the sample under study. 



C. Center-of-mass level representation of the 
polymer chain 

In the soft colloidal particle representation the projec- 
tion operator targets the center-of-mass of the polymer. 
The field variable {n = I, N = I, a ~ cm) is simply 
defined as 



g (R, P, t) = [6 (r„„(i) -R)5 (p„„(<) 



(19) 
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Applying the projection operator in the new field variable 
to the Liouville equation, where U = and Wij ^ 0, 
leads to the generalized Langevin equation 

^g(R,P,i) = - J ds J dR' J dP'M(R,P,R',P') 
x.g(R',P',(t-s)) +^^(R,P,t) ,(20) 

which reduces, following the procedure briefly outlined 
in Section [Tir^ to 

+ Fi.(f). (21) 

In the overdamped regime, 

W^^=FiW, (22) 

where Csoft is the friction coefficient for the col- 
loidal particle, (soft = /3/3 dr(F,„(i) • Ff^{t - 
t)). Eq.(|22| obeys the fiuctuation-dissipation relation 

{F,m{t) ■ F7m{t')) = St-t'6kBTCsoft- 

The choice of the field variables in the projection op- 
erator defines the length scale of coarse-graining and the 
variables in which the resulting Generalized Langevin 
Equation is expressed. Because the derivation just pre- 
sented depends on the basic assumption that the correla- 
tion function of the bath variables are short lived in the 
presence of heavy particles, and correction terms repre- 
sented by the memory functions are minimized when a 
clear separation of timescale is observed between the slow 
variables in the projection operator and the fast variables 
that are averaged out, this criteria provides a way of se- 
lecting the relevant lengthscales for the coarse-graining, 
when dynamical properties are under study. 

For example, as far as polymer dynamics is concerned, 
we know that for times longer than the longest Rouse 
correlation time, tj^ w R^/D, polymer internal dynam- 
ics is fully relaxed and the monomer dynamics follows the 
motion of the center-of-mass, which is long lived. This 
suggests that the center-of-mass coordinates are a good 
choice to represent the projected slow dynamics for time 
t » tr. This reasoning holds for both unentangled 
and entangled polymer dynamics as the longest relax- 
ation time, after which free diffusion and Brownian mo- 
tion set in, is tr with the proper diffusion coefficient, i.e. 
for unentangled chains Dunent oc and for entangled 
chains -Dent oc 7V~^. 



should predict identical dynamics for the polymer center- 
of-mass, shows that the two equations differ because of 
the presence of the intramolecular free energy in the 
monomer description, which is absent in the soft colloidal 
approximation, and because of the different friction co- 
efficients in the two representations. 



A. Free energy rescaling 

The elimination of degrees of freedom increases the en- 
tropy of the system, as every coarse-grained state corre- 
sponds to a number of preaveraged microstates. In an ex- 
treme picture we can imagine that the preaveraging due 
to the coarse-graining procedure is in effect transform- 
ing the energy of the system, expressed for example in 
the Liouville equation by an Hamiltonian, into a free en- 
ergy in the corresponding Langevin equation. While the 
Hamiltonian contains kinetics and potential energy, the 
free energy includes an entropic contribution due to the 
preaveraged microstates for each coarse-grained state. 

As far as the free energy correction is concerned, the 
system described by the larger cutoff lenghtscale is the 
one where the level of coarse-graining is most extensive 
and the highest entropic correction has to be included. 
This correction can be calculated from the comparison 
of the two equations. Because the system described at 
the monomer level is exploring in time the intramolec- 
ular energy states of the configurational landscape, its 
dynamics is slowed down with respect to the colloid rep- 
resentation where intramolecular degrees of freedom are 
not present. To take this effect into account we calculate 
the correction that has to be included in the soft-colloid 
representation to take into account the time spent by the 
atomic system to explore the internal degrees of freedom. 

Consistent with the monomer-level model adopted in 
our study and with UA MD simulations, the polymer 
is described as a collections of beads, or friction points, 
connected by harmonic springs. Each bead corresponds 
to a CHx moiety, with a; = 2 or 3, depending if the unit 
is imbedded in the chain or is terminal. This model has 
been shown to provide a realistic representation of the 
dynamics of numerous polymeric systems with different 
chemical structure. [571 [551 liO l [5 ^ 1^5] 

The intramolecular potential is defined as 

f^W = ^E ' (23) 



IV. 



ANALYTICAL RESCALING OF THE 
COARSE-GRAINED DYNAMICS 



The two Langevin equations, Eqs.(15) and (22), dis- 



play the two levels of coarse-graining of the macromolec- 
ular liquid, which are adopted in this paper. The compar- 
ison of the two equations, which in the longtime regime 



with U{r) not to be confused with the semifiexibility ma- 
trix of Eq.(18). Here A is the connectivity matrix of 
Eq.(ll6|), which represents the structure and local flexi- 



bility of the polymer [66l ISl] , the position of bead i in 
a chain of TV beads or united atoms, and li = r^+i — 
the bond vector connecting two adjacent beads. 

The statistically averaged internal energy for one 
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molecule consisting of N monomers is given by 



U 



= N / J/e"^ 



Ar 



dr — 



37V 



r Are 



rr" Ar 



(24) 

After solving the integral by normal mode transforma- 
tion, as reported in Appendix I, this model predicts the 
average energy dissipated in the internal modes to be 
{U/ikBT)) = 37V/2. The soft-colloidal representation, 
instead, has no internal degrees of freedom. 

The simulation time t, as measured in the MS simula- 
tion of the coarse-grained system, translates into the real 
time t after including the rescaling due to the energy, 
which is reduced by the amount of energy dissipated in 
the fluctuations due to internal degrees of freedom. [55] 
For our model 



, m 3 



kBT2 



(25) 



with the particle mass, m, and size Rg. This rescaling 
slows down the coarse-grained dynamics, but only par- 
tially accounts for the observed phenomenon because the 
rescaling of the friction needs to be included. 



is Csoft/ {N(m), which yields the rescaled mean-square 
displacement 



{AR\t))=6D 



MS Csoft 



(27) 



The thermodynamic conditions of the system under 
study, i.e. density and temperature, and its molecular 
structure, i.e. the radius-of-gyration, enter the equation 
above both directly through the definitions of the friction 



coefficients, Eqs.(44| and (49) and indirectly through the 



mesoscale simulation from which the diffusion coefhcient, 
D^^^ , is measured. 



To solve Eq.(27) we start from the definition of the 
monomer friction coefficient, Cm, which is given in the 
Markov limit by the memory function 



1 

N 



N 

E 

a, 6=1 



(28) 



^a.b{i) is the function that describes the correlation, 
through the polymer chain between monomers a and 6, 
of the random forces generated from the random colli- 
sions of the surrounding molecules undergoing Brownian 
motion, with Ii58n60i 



B. Monomer friction coefficient 

The rescaling of the friction coefficient is calculated 
considering the friction of the polymer center-of-mass in 
the monomer/UA representation, and comparing the re- 
sult with the friction of the cm of a soft colloidal particle. 
The expression for each of the friction coefficients is de- 
rived from its definition as the integral of the memory 
function contribution to the Generalized Langevin Equa- 
tion (GLE) in the two levels of representation. 

The effect of coarse-graining the Liouville equation, 
or projection onto the slow degree of freedoms, is the 
appearance in the Langevin equation of the dissipation 
terms, given by the random force and the friction coef- 
ficient. Systems with different levels of coarse graining 
have different friction and, as a consequence, different 
diffusion coefficients. 

For a particle in a liquid, the center of mass mean- 
square displacement is defined as 



(Ai?2(t)) ^&Dt 



(26) 



with D the diffusion coefficient. For a polymer, the cm 
diffusion coefficient is given by Z? = ksT /{NC,rn), where 
Cm is the friction coefficient of a monomer, while for a 
liquid of soft colloidal particles Dsoft — ksT / C,soft, with 
Csoft the friction coefficient of the colloidal particle. The 
two should be identical in the long-time limit, but they 
are not, as the diffusion coefficient obtained from MS MD 
simulation is much larger (much faster dynamics) than 
the one obtained from UA MD. The correction factor 
to scale down the MS MD diffusion coefficient, D^^^, 



ra,b(t) = ^pjdrj dv'g{r)g{r')F{T)F{r')v-v' 

X j dKS%{R-t)S'^{\Y~v' + n\-t), (29) 

where g{r) = h{r) + 1 is the monomer radial distribution 
function, F{r) is the total force exerted by all the matrix 
polymer on the monomer, and S^{r\t) is the projected 
dynamic structure factor of the matrix fluid surround- 
ing the polymer. The unit vectors f and f' define the 
directions of the total exerted forces. The derivation of 



Eq. ( 29 ) is not completely new and is briefly reported in 



Appendix II. 



Eq. ( 29 1 rests on the approximations that the fluid is 



isotropic and that many-body correlation functions can 
be described with good accuracy as products of pair 
distribution functions. The solution of this equation is 
sometimes carried on by introducing a mode-coupling 
approximation. [55^61] however we follow a different pro- 
cedure. The dynamic structure factor, which is ruled by 
the projected dynamics, is approximated by its real dy- 
namics counterpart, S^{r;t) w S{r;t), simplifying the 
solution of Eq.(29). This is an acceptable approximation 



when the Langevin equation is expressed as a function 
of the slow variables [65 and holds for our system in the 
long-time, diffusive regimej62]. 

In order to separate the spatial coordinates of S{\r — 
r' -I- R|;t) in Eq. (29) it is convenient to use the Fourier 
transform 



Sir;t) 



1 



(27r)3 



e'^'^Sik; t) dk 



(30) 
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where the dynamic structure factor is calculated in re- 
ciprocal space as the sum of intra- and inter-molecular 
contributions 



07 



+ (31) 



Here uJaf{k,t) is the time dependent intramolecular 
probability distribution functions for monomers a and 7, 
on the same molecule, to be separated by a reciprocal dis- 
tance k, while haj(k,t) is the corresponding intermolec- 
ular contribution. 

Given that the dynamics on the global scale is driven 
by the polymer diffusion, the intramolecular probabil- 
ity distribution function in reciprocal space can be ex- 
pressed, in the limit of large lengthscales, k < 1/Rg, as 



UJa-yik; t) « CXp 



71 



6 



exp[-A:2£)t] ^ (32) 



where D is the polymer center-of-mass diffusion coef- 
ficient and I = ^^^X^ili l^il ^^^^ average segmental 
length. 



Because in Eqs.(28) and (29) the order of the sum- 



mation and time integrals can be changed, the double 
summation reduces the inter- and intramolecular dis- 
tributions to their averages over the bead distribution. 
The site-averaged intramolecular probability distribu- 
tion function, u}o{k)^ is well approximated by the Debye 
formula, EH 



Wo 



1 ^ 27V(e'='-"' + k^Rl - 1) 



or by its Pade' approximant 



u!o{k) 



N 



1 + /fc2^2 



(34) 



The site-averaged intermolecular probability distribu- 
tion is defined by the Ornstein-Zernike equation 



^Uk)c{k) 



pc{k)ujo{k) ' 



(35) 



where c{k) is the direct correlation function. At the 
monomer level we follow Curro and Schweizer's PRISM 
thread approach, [53l [54] where the polymer chain is mod- 
eled as a thread of vanishing thickness, c(fc) ss cq, with 
Co = -(1 - 2ep/Rl)/{2N^Pckep/Rl)- Substitution of co 



and Eq. ( 34 ) into Eq. iSbh gives 



hik) 



ho 



(l + fc2e2)(l + fc2e2) ' 



(36) 



where ho = h{k = 0) = (Cp/Cc ~ ^)/Pch is related to the 
compressibility of the system. [211 [13] 



Because in the large length scale regime, of interest 
here, the relaxation of the liquid is dominated by the 
polymer diffusion, the dynamic structure factor is ap- 
proximated as 



S{k;t) « S'(fc)exp [-k^Dt] 



(37) 



Finally, after introducing the integral representation of 
the delta function 



J dR e'^'^^'+^'^ = (2^)3 5(ki + ka) 



(38) 



the last integral in Eq.(29) simplifies to 

J dILS{R;t)S{\v - r' + R\;t) = J dki J dk2 

S{ki;t) S{k2;ty^^'^'-''^ J dR e'^^^'+^'^ . (39) 

Because the functions wo(fc) s-nd h{k) are even with re- 
spect to fc, the equation reduces to three contributions: 
the first is due to intramolecular interactions iUo{k), the 
second includes the cross product ujo{k)h{k), and the last 
is due to the intermolecular contribution h^(k). This 
leads to the following expression 

I dRS{R;t)S{\r - r' + R\;t) = J dk 



(wo(fc) -I- 2ph{k)ujo{k) -t- p'^h^(k))e 



-2k^Dt i k(r-r') 



(40) 



Because we are assuming that monomers are interact- 
ing through a hard core potential, which is consistent 
with the PRISM thread model, [ill [Si] the force is a delta 
function and therefore 



g{r)F{r)^g{d)r^5{r-d). 



(41) 



where d is a hard core diameter, identical for any CH2 
bead in the chain, in the spirit of the UA-MD description 
and PRISM approach. When we compare our equations 
with data of experimental or simulated systems, where 
monomers interact through a Lennard- Jones potential, 
the latter has to be mapped onto a hard-core potential 
with the effective diameter, d.[55] 

The final expression for the monomer friction coefh- 
cient is given by 



■pg\d) {pDY 



Cm - ^ 

x{j[u:o{k),uJo{k)]+2pJ[ujo{k),h{k)]+ p'J[h{k)Mk)]] 



(42) 



with the function 
J[a(fc),/3(fc)] = I dr I dr' 



',,^Hk\r-^,,, 

fc|r-r'| 
X J(r - d)5{r' - d) a{k)f3{k) . (43) 
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The solution of Eqs.(42) and (43 1 is given by a lengthy 
but analytical expression, which is a function of the 
molecular parameters, ^p, i?„, thermodynamic parame- 



ters, p, /?, the diffusion coefficient, _D, and of the hard- 
core diameter d, as 



15\/2-F40— + 12^2 ^ 



Ra 



Rr. 



pirNhQ 



3V2{Rj - 2^2)2 



a 

4o3li .1 ^pY\ , 4V2d3i?4 f 5 _ 14 f ^lV + 2 f 4^ V I + 3d2i?5 ( 5 - 14 ^ ' ' 



12^2^; 



(44) 



42d^i?^(l-2(^ 



4^2 



42^26 ^'■Cl 1 



^'"'^°12(i?2 2.2)3 



AOd^Rl + l5V2d^Rl ~ 24:V2d'^Rl£.l - 144d='i?^e, 



c2 
"gSp 



+6V2d^Rl (2 - 9 (^l)'"! + 12i?2eJ U 



9 -8^,^4 
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This expression is general and holds for any homopoly- 
mer melt represented as a collection of identical beads in- 
teraction through a hard-core potential of range d. The 
value of d is specific of the monomeric structure of the 
homopolymer. 



C. Friction coefficient for a liquid of interacting 
soft colloidal particles 

The friction coefhcient for a point particle interacting 
through a soft repulsive potential is much smaller than 
the friction of the macromolecule before coarse-graining. 
In fact, the friction coefficient of an object can be esti- 
mated using Stokes' formula where ^ = GirrirH, with i] 
the fluid viscosity and rj^ the hydrodynamic radius. The 
latter can be evaluated from the surface area of the ob- 
ject exposed to the solvent, which can be estimated by 
"rolling" a solvent molecule on the object. It is evident 
that the surface available to the solvent in a bead-spring 
representation of a polymer is much higher than the sur- 
face available to the solvent for a point particle interact- 
ing through a soft, long-ranged potential. 

To calculate the friction coefhcient for a soft colloidal 
particle, we start from the Generalized Langevin Equa- 
tion that describes the time evolution for the position 
coordinate of the molecular center-of-mass, i.e. Eq.(22) 



where the friction coefficient for soft particles is given by 



t,soft 



(/3/3) Pch / dt dr dr'g{r)g{r')F{r)F{r') 



xr • r 



'J dRS{R-t)S{\r-r' + 



t). 



(45) 



form of the density prefactor. In reality the form of the 
pair-distribution function, g{r), the force exerted by the 
surrounding molecules on the tagged chain, F(r), and the 
dynamic structure factors, 5(R, t), are different quanti- 
ties in the monomer and soft-colloid representations. 

We assume that the dynamic structure factor in recip- 
rocal space has the form 



S{k; t) « S{k) e-"'" = (1 + Pckh'^'ik)) e 



(46) 



where h'^'^{k) is the center of mass total pair correlation 
function 



h'^'ik) = ho 



1 + k^Rl/2 

1 + k^ep 



(47) 



with ho = ^ ^)/Pch, as defined in the previous 

section. Eq.(|47| is just the Fourier transform of Eq.([I]). 



Eq.(46) indicates that in the long-time regime, which is of 



Eqs.(29) and (W5| look identical, with just a different 



interest here, the relaxation of the liquid is largely driven 
by the center-of-mass diffusion, while internal dynamics 
and local modes of motion are already fully relaxed. This 
is a reasonable assumption given that the lengthscale of 
our treatment is the overall polymer dimension, and no 
structural or dynamical information is retained on the 
local scale. 

To perform our calculation we need to defined an ap- 
proximate analytical form of the effective force. To do 
so we adopt the simplified form of h'^^{r), Eq.([2]). Then, 
we reduce further the expression by neglecting the small 
attractive component of the potential. Finally we ap- 
proximated the real potential, v{r), with its mean- force 
counterpart w(r) « — fc^Tln [h{r) + 1] properly rescaled. 
The real potential, calculated through the HNC approx- 
imation as described in the following section, is a com- 
plicated function of h{r). However it can be related, in 
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an approximated way, to the simpler potential of mean 
force through the equation 



v{r) 



w{0) 



w{r) 



(48) 



where v{r) « ^/3w{r) for all the samples considered in 
this study. These approximations define the force F{r) 
and the pair distribution function, g{r) = h{r) + 1, en- 
tering the equation for the friction coefficient. 

The resulting expression for the friction coefficient of 
the soft colloidal particle is expressed as a function of the 
diffusion coefficient D, /3, p^h and the two length scales 
Rg and as 




3 1183 , 679\/3 ^ ,2 



(49) 



This expression is an approximated analytical form for 
the friction coefhcient of a soft colloidal particle. 



V. MESOSCALE SIMULATIONS 

Here we present numerical calculations to illustrate 
and discuss the rescaling procedure of the preceding sec- 
tions. We first perform MS MD simulations of the coarse- 
grained polymer liquid, where each chain is represented 
as a soft colloidal particle, centered at the center-of-mass 
of a chain, and interacting with the surrounding parti- 
cles through a soft repulsive potential of the order of few 
times the chain dimension, Rg. The simulations of the 
soft colloidal liquid produce dynamical properties that 
are accelerated due to the soft nature of the potential 
in the coarse-grained representation. These properties 
are rescaled following our procedure, and then compared 
with existing data, when they are available. 

In a previous paper we briefly presented calculations 
of the rescaled dynamics for a variety of systems in- 
cluding UA MD simulations and experimental data of 
PE diffusion available in the literature, that we use 
to test the accuracy of our procedure. We selected 
UA MD simulations as our test (see Table |l| because 
they have been shown to reproduce with a high level 
of accuracy the dynamical properties of PE melts, such 
as diffusion and viscosity. PTTI, We also compared 
predictions of rescaled MS-MD simulations with ex- 
periments for samples with temperature T = 509 K, 
monomer site density p — 0.0315302 sites/A'^ and N = 
36, 72, 106, 130, 143, 192, and 242. 39, 48 51,. Our 
MS-MD simulations, properly rescaled, provided good 
quantitative predictions of the diffusion coefficient for 
those systems. |3lj 

In this paper we use those same systems to illustrate 
our procedure. Moreover we present new results for PE 



samples, not present in the literature, to underline the 
predictive power of the theory, where no calibration curve 
is necessary. Once a system is selected, its structural and 
thermodynamic parameters are defined and are used as 
input to the MS MD simulation so that the whole proce- 
dure is free of adjustable parameters, with the exception 
of the parameter d that is fixed for PE once and for all 
samples. [14, 45-47 

Systems that we simulated include liquids of chains 
with increasing degree of polymerization, as described 
above. As the molecular weight of the polymer increases, 
the systems cross the threshold from unentangled to en- 
tangled dynamics. For entangled systems the dynamical 
rescaling approach that we propose is modified to include 
a one-loop perturbation that accounts for the presence 
of entanglements. Simulations of soft colloidal liquids 
are performed for entangled systems, and the rescaling 
applied to predict diffusion. 



TABLE I: Polyolefin melts UA MD simulation parameters 



System 


T[K] p[sites/A^] 






PE 30" 


400 


0.0317094 


63.5695 


81.7544 


PE 44" 


400 


0.0323951 


110.3197 


127.6856 


PE 48'' 


450 


0.0314487 


111.0832 


140.8119 


PE 66" 


448 


0.0328993 


177.5348 


199.8812 


PE 78'' 


450 


0.0321465 


205.9221 


239.2607 


PE 96" 


448 


0.0328194 


281.7989 


298.3301 


PE 122'' 


450 


0.0325479 


346.2655 


383.6526 


PE 142'' 


450 


0.0326600 


420.7070 


449.2852 


PE 174'' 


450 


0.0327680 


525.1816 


554.2975 


PE 224'' 


450 


0.0328835 


690.5038 


718.3791 


PE 270'' 


450 


0.0329520 


856.4648 


869.3342 


PE 320'' 


450 


0.0330034 


980.1088 


1033.4158 


" from Refs. 


45ti47]: '' from Ref. [HITS] 





Details about our MS MD simulations have been re- 
ported in previous papers of ours and will not be re- 
peated here. UHl [301 131] Briefly, MS MD simulations were 
implemented in the microcanonical (N, V, E) ensemble 
on a cubic box with periodic boundary conditions. We 
used reduced units such that all the units of length were 
scaled by Rg (f = r/Rg) and energies were scaled by 
ksT . Temperature and radius-of-gyration were utilized 
for dimensionalizing the results obtained from the MS 
MD simulations, after they were performed. 

The number of particles, i.e. polymer chains, in our 
simulations varies from 1728 [N = 40) to 85184 {N = 
1000) depending on the system. This number is deter- 
mined by the box size, which is larger than twice the 
range of the potential, and by the liquid density. The 
potential is long-ranged, due to the many-body effects 
entering through the OZ equation. 

Each simulation evolves for about 50, 000 computa- 
tional steps. For the entangled melts the potential is 
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longer ranged than for the unentangled systems, and 
therefore it is cut at larger distances requiring a big- 
ger box size. The reduced density used in simulations, 
pstm _ pj^z where p is the site density, varies around 
1 for unentangled melts, and exceeds 2 for weakly en- 
tangled melts. A typical MS simulation takes between 2 
hours [N = 40) to 4 days [N = 200) on one CPU work- 
station, while using the code that works in parallel, the 
computational time is further reduced. 



A. Interparticle potential 

The pair potential acting between two effective coarse- 
grained units is formally derived from the colloidal rep- 
resentation of the liquid, specifically h{r), using an 
hyper-netted-chain (HNC) closure approximation to the 
Ornstein-Zcrnike equation. |56l This approximation is 
known to work well for liquids of particles interacting 
through a soft potential. [69| The potential input to the 
MS MD simulation, w'^'^(r), is derived from the total cor- 
relation function for the soft colloidal representation of 
the liquid, h'^'^[r), defined in the limit of long chains, 
N ^ CO, as in Eq.([2|. The potential is calculated using 
the hypernetted chain approximation as 

Pv'^'ir) = h^'ir) - ln[l + /i""(r)] - c'"'{r) . (50) 

Here the direct correlation function, c^^(r) is given in 
reciprocal space in terms of h'^'^{k) as: 



c-(fc) = 



l + Pchh'^'ik) ■ 



(51) 



It is important to define the correct potential acting 
between the coarse-grained units to achieve a realistic 
representation of the large scale properties of a system 
through MS MD. Because coarse-grained potentials re- 
sult from the mapping of many-body interactions into 
pair interactions, through the averaging over microscopic 
degrees of freedom, they are parameter dependent. Dur- 
ing the coarse-graining procedure, the potential acting 
between microscopic units, which is given by the Hamil- 
tonian of the system, reduces to an effective potential, 
which is a free energy in the reference system of the mi- 
croscopic coordinates. The coarse-grained potential so 
obtained contains contributions of entropic origin due to 
the microscopic, averaged-out degrees of freedom and is 
therefore state-dependent. This can be observed in the 
form of the total correlation function between coarse- 
grained sites, Eq.([2]) from which the potential is derived. 
The correlation function explicitly includes the structural 
and thermodynamic parameters of the polymer, i.e. the 
radius-of-gyration, density screening length, and number 
density. The temperature enters directly through Eq.( 50 ) 



B. Results from mesoscale molecular dynamics 
simulations 



Before entering the details of applying our rescaling 
approach we focus on the "raw" dynamics obtained di- 
rectly from the MS MD simulations. Fig. [3] displays the 
mean-square displacement for the MS MD simulation of a 
polyethylene melt with N = 44. At short times the iner- 
tial term in the Langevin equation is dominant as the par- 
ticles undergo ballistic dynamics, while in the long-time 
regime the system crosses over to diffusive dynamics. The 
diffusion coefficient is higher than the value measured 
in UA MD simulations, and the transition from ballistic 
to diffusive regime happens after about 5000 simulation 
steps (dot-dot vertical line on top panel of Fig.js]), which 
corresponds to a distance of roughly 30i?g. Such a large 
distance reflects the fact that in MS MD simulation the 
point particles interact through a very soft potential and 
the density is also very low. Because the particle has to 
" collide" many times to undergo the crossover to diffusive 
dynamics the latter takes place at a large lengthscale. 

Moreover, the bottom panel in Fig. [3] displays the 
velocity correlation function {{v{t) — w(G))^) and shows 
that, consistent with the mean-square-displacement, once 
more the inertial term becomes negligible at the same 
crossover time that the diffusive regime sets in. Since 
our MS MD simulation are performed at equilibrium, 
the avarage kinetic energy per particle (mf^(i)/2) = 
3/2 ksT, and therefore 

lhn{iv{t)-viO))') = lirn [2{v'it))~2{vit)v{0))] 



QksT /m . 



(52) 



and indirectly through the molecular parameters, such as 



The dynamical transition is displayed as a dashed line 
on the bottom panel of Fig. [3) taking into account that 
our simulations are in reduced units and m = 1, kgT ~ 
1. The figure shows that the velocity autocorrelation 
function reaches it's asymptotic value at about the same 
time as the diffusive regime sets in for the mean-square- 
displacement. 



VI. APPLICATION OF THE RESCALING 
PROCEDURE 



As stated above, the accelerated dynamics that is a 
consequence of the coarse graining of the system can be 
rescaled by taking into account the two main effects of 
the procedure, namely the change in the entropic contri- 
bution to the free energy of the system due to the aver- 
aging of the internal degrees of freedom and the change 
in the friction coefRcient due to the different shapes of 
the molecule in the two different levels of coarse-graining. 
The difference in shape relates to the change in the molec- 
ular surface available to the surrounding molecules, and 
to the correlation of the random forces generated by in- 
termolecular collisions. 
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FIG. 3: Top panel: cm mean-square displacement (solid line) 
from MS MD simulation in reduced units as a function of sim- 
ulation time steps for the PEAA melt sample. The slopes for 
the ballistic and the diffusive regimes are shown as dashed and 
dot-dashed lines correspondingly. Bottom panel: cm veloc- 
ity time correlation function showing when the inertial term 
becomes negligible. The asymptotic value of 6 kBT/m is de- 
picted as a dashed line. 



The first rescaling is given by the inclusion a posteriori 
of the internal degrees of freedom, averaged out during 
the coarse-graining procedure, as a correction term in the 
free energy of the system, which accounts for the differ- 
ence in entropy. The energy correction affects the time 
of the measured dynamics as the change from a bead- 
spring description to a soft-colloid representation leads 
to the rescaling of the time reported in Eq.( 25 ), also tak- 



ing into account the fact that because the potential is 
expressed in normalized quantities, the simulation runs 
using reduced units of energy, fc^T = 1 and the normal- 
ized length r/Rg. 

The second rescaling of the dynamics is calculated 
starting from the ratio between the friction coefhcients 
in the two coarse-grained representations, as described 
in Eq.(l27l). 



A. Calculations of the monomer friction 
coefficient, 

Because our formalism maps the Lennard- Jones liquid 
described by the UA MD simulation into a liquid of poly- 
mers interacting through a monomer hard-core repulsive 
potential, it is necessary to define an effective hard-core 
diameter, d. This is done by requiring the friction of the 
chain with iV = 44 to follow the expected scaling be- 
havior for the diffusion of an unentangled polymer chain. 



D = k bT / [N C,rn) ■ Since all except two of the atoms 
in our PE chains are CH2 monomers, we assume that 
the potential is identical for all the units along the ho- 
mopolymer chain. Moreover we assume that the range 
of the repulsive interaction, d, is independent of liquid 
density. 

Among the different samples, we selected the chain 
with = 44 to optimize d, because this sample fol- 
lows unentangled dynamics while the polymer is long 
enough to obey the Gaussian intramolecular distribu- 
tion of monomer positions, which justifies the analytical 
form of the intramolecular structure factors used in our 
formalism. FigJH] displays the monomer friction coeffi- 
cient, from Eq.(44), expressed as the dimensionless quan- 
tity Df3(^„i, as a function of hard sphere diameter d, for 
polyethylene melts of three different degrees of polymer- 
ization. The 1/A^ scaling is reported as a dot-dashed line 
in the figure. 

In these calculations, the numerical values of A^, p and 
Rg were taken from the data of the UA MD simulation 
against which the proposed approach is tested. The value 
of the radial distribution function at the contact was set 
to g{d) = 1/2, which is the conventional value assumed in 
the PRISM thread theory for polyethylene chains. This 
value is intermediate between zero and the first solva- 
tion shell value. The optimized hard-core diameter for 
A^ = 44 is d = 2.1A, which is an intermediate value 
between the bond length, I = 1.54 A, and the Lennard- 
Jones cr-parameter, a = 3.95 A, in the intermonomer 
potential of the UA MD.gH g?! The unentangled scal- 
ing is fulfilled for PE30 at d = 2.07 A and for PE96 at 
d = 1.96 A, which are close to the one for PE44. Because 
the PE96 sample has a degree of polymerization that is 
close to the entanglement value of A^e = 130, its dynamics 
is likely to be in the crossover regime where the effect of 
entanglements start to be felt, since the transition from 
the unentagled to the entangled dynamics is very broad. 

Table |lT] displays the numerical values of the dimen- 
sionless monomer friction coefficient, DfiQm, for poly- 
meric liquids with different degree of polymerization, 
A^, across the untentengied-to-entangled transition. For 
both unentangled and entangled systems the hard sphere 
diameter has been fixed to the value of the unentangled 
ones, d — 2.1A, so that the intermolecular monomer po- 
tential is not changed as a function of A^. While for un- 
entangled systems the monomer friction coefficient was 
calculated from Eq.(44), for entangled chains we adopted 



a perturbative approach to account for the effect of en- 
tanglements. Let's denote 



DPCm^Jip.N.Rg^d), 



(53) 



where for unentangled systems J{p,N,Rg,d) ~ N~^. 
Following a one-loop perturbation, and including the def- 
inition of the diffusion coefficient for a macromolecule 
comprised of A^ monomers with the monomer fric- 
tion coefficient D = {f3(jnN)~^ , the normalized and per- 
turbed friction coefficient becomes 



D13C = NDPQmJip, N, Rg, d) = N{Df3Uy^ 



(54) 
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FIG. 4: Plot of DfiC,m as a function of hard core diameter d. 
From left to right, the three panels show curves for polyethy- 
lene melt with iV = 30, iV = 44, and iV = 96. 



The one loop perturbation is in the spirit of the repta- 
tion model where both the chain reptating and the chains 
involved in the entanglements relax with the same diffu- 
sive mechanism: each brings a scaling contribution, 
which is the trademark of polymer Brownian motion. In- 
terestingly, in our model the diffusion coefficient of entan- 
gled polymers under certain fixed monomer density and 
temperature shows apparent scaling exponents different 
from —2 (see Fig. 6|. The resulting scaling exponents 
emerge cumulatively from output of mesoscale simula- 
tions and both steps of rescaling. 



Because Eq. ( 54 ) applies only when the systems are en- 
tangled, to predict the diffusive behavior of new samples 
it is necessary to estimate a priori the crossover degree 
of polymerization, N,,. Several methods have been pre- 
sented in the literature to estimate from thermody- 
namic conditions and molecular parameters. j73] Those 
methods provide similar values of A^e- Moreover, the 
expressions for the unentangled and entangled frictions, 
Eqs.(44| and (54), predict values that differ only slightly 
in the crossover region, as it is shown in Fig. |6] in this 
paper. In this way, selecting the unentangled expres- 
sion to represent entangled systems, or viceversa, in the 
crossover region would result in small inconsistencies in 
the calculated diffusion coefficients. 

Table In] includes intra- (-D/3C„i(wQ-y)) and inter- 
molecular contributions to the monomer 
friction coefffcient, as well as the self intramolecular con- 
tribution, D/3(^rn{^aa)- In general, the calculated to- 
tal friction is comparable in magnitude to the self in- 
tramolecular contribution, DP(^mi^aa)- Moreover, the 
total intramolecular contribution, D/3(^m(j-^a,-y), is of the 
same order of magnitude of the itermolecular contribu- 
tion, D(3(^miha,f), but with the opposite sign, which is 
reasonable as the liquid is almost incompressible. This 
result shows that the conventional approximation of re- 



placing the structure factor, Sa,j{k), by the single chain 
analog, uJa,'y{k), can lead to errors in the evaluation of 
the memory function for macromolecular liquids. [71| The 
table also displays the value of the dimensionless fric- 
tion coefficient ND/3C,n, which for unentangled systems 
should be « 1. As expected, we see deviation from the 
unentangled behavior in the very short chains and in the 
crossover to entangled dynamics at iV sa 100. 



TABLE II: Monomer friction coefficient contributions with 
hard sphere potential for d — 2.1 A 



System 






) D/3Cm{hc-y) 






PE 30 


0.03378 


0.08484 


-0.04883 


0.03601 


1.0804 


PE 44 


0.02744 


0.06244 


-0.03991 


0.02253 


1 


PE 48 


0.03101 


0.07984 


-0.04904 


0.03080 


1.4782 


PE 66 


0.02539 


0.05703 


-0.03823 


0.01880 


1.2409 


PE 78 


0.02632 


0.06287 


-0.04179 


0.02107 


1.6439 


PE 96 


0.02239 


0.04766 


-0.03325 


0.01441 


1.3834 


PE 122 


0.02400 


0.05533 


-0.03831 


0.01701 


2.0757 


PE 142 


0.02245 


0.04967 


-0.03501 


0.01466 


2.0815 


PE 174 


0.02193 


0.04827 


-0.03438 


0.01389 


2.4170 


PE 224 


0.02133 


0.04664 


-0.03359 


0.01304 


2.9219 


PE 270 


0.02039 


0.04344 


-0.03164 


0.01180 


3.1849 


PE 320 


0.02184 


0.04956 


-0.03585 


0.01371 


4.3874 



B. 



Calculation of the friction coefficient of a soft 
colloid, (soft 



Starting from Eq. ( 49 1 we calculated the friction coeffi- 
cient, DPQsofti for polymer liquids represented as soft 
colloidal particles. Table |III| shows the dimensionless 
friction coefffcient for several systems. The molecular 
parameters, TV and -Rg, and the thermodynamic condi- 
tions of density p and temperature T, are taken from the 
UA MD simulations, see Table |Tj The dimentionless fric- 
tion coefficient for these systems is Z?/3Csoft(l + h'^'^) ~ 
0.002 - 0.01, while we would expect DjiCsoitll + h"") « 1 
for unentangled systems, see Table [111} These data show 
that the theoretically calculated friction coefficient (with- 
out rescaling) for the soft colloidal systems greatly un- 
derestimate the friction coefficient, as also observed in 
the MS MD simulations, and hence give rise to accel- 
erated dynamics as discussed previously. It also shows 
that intra- and inter-molecular contributions to the fric- 
tion coefficient are comparable in magnitude; both of 
them need to be taken into account when calculating dy- 
namical properties of polymer melts. 
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TABLE III: Soft colloids friction coefficient contributions 



System 








rrj oU 


n C\A A OVQ 
U.U44Z / 


-U-Uzyyoz 


U.Ui4o4i 


rill 44 


U.U2U ( oy 


-U.Uiz441 


U.UUoozo 


Jr^rL/ 4o 


U.Uz4Doy 


n ni p;oQO 
-U.Ulozoy 


u.uuyo4y 


Jr^Jl/ DO 


u.uizDiy 


-U.UUDDoy 


u.uuoyou 


r^il/ / 


u.uizuiy 


-U.UUDZ04 


U.UUo / DO 


PE 96 


U.UUOUtJtJ 


U.UUO I 1^ 




PE 122 


0.007423 


-0.003334 


0.004089 


PE 142 


0.006159 


-0.002612 


0.003546 


PP 17zL 


0.005218 


-0.002108 


0.003109 


PE 224 


0.004298 


-0.001648 


0.002650 


PE 270 


0.003660 


-0.001351 


0.002310 


PE 320 


0.003521 


-0.001288 


0.002233 


C. 


Results from the rescaling procedure. 



Comparison with simulation and experimental data 

Some of the results reported in this section were al- 
ready briefly presented in our short paper. [3T] Our dis- 
cussion here makes use of some of those data as a starting 
point to illustrate with an example the details of the pro- 
posed rescaling procedure and highlight its strengths and 
weaknesses. 

In order to rescale the unrealistic fast dynamics of MS 
MD simulations we applied our rescaling procedure and 
compare the predicted dynamics with data from UA MD 
simulations and experiments. We use as input param- 
eters the thermodynamic conditions and molecular pa- 
rameters of each sampl e un der study. The rescaling pro- 
cedure is given by Eq.(27), where is the diffusion 



coefficient from the MS-MD simulation, the soft colloid 



friction coefficient is calculated using Eq.(49), and the 
monomer friction coefficient is given by Eq. ( 44 ) for un- 



entangled chains, and Eq.(54|with Eq.(44) for entangled 
ones. Eq.(27) depends on the temperature and density 
of the system investigated and on its molecular radius- 
of-gyration. 

Indirectly those parameter enter our procedure 
through the diffusion coefficient from mesoscale simula- 
tions, I?*^"^. Specifically temperature enters through the 
rescaling of the time, as the time step in the mesoscale 
simulation is adimensional and becomes dimensional once 
it is rescaled by the energy, following a well-established 
procedure. Moreover, thermodynamic and molecular 
parameters enter indirectly through the soft potential, 
Eq. ( 50 1 , which is parametric and includes density, tem- 



perature, and the molecular radius-of-gyration. 

Finally, thermodynamic parameters enter directly 
through the definitions of the friction coefficients in the 
monomer and soft-sphere descriptions, Eqs.(44) and (49) 
respectively. Specifically, the monomer friction coeffi- 
cient is a function of p, iV, Rg, plus a hard sphere di- 



ameter, d, which is used to map the Lennard- Jones po- 
tential of the united atom simulation onto a repulsive 
hard-core potential with an effective bead diameter. The 
hard sphere diameter d is assumed to be independent of 
the thermodynamic conditions, for the range of temper- 
ature and density simulated here, and constant for all 
the monomers in the homopolymer chain. The criteria 
of choosing numerical value for d have been already ex- 
plained and discussed. 

In an analogous way, the soft-sphere friction coeffi- 
cient depends on the chain number density, which re- 
lates to the monomer number density through N as 
Pch = p/N , and Rg is the radius of gyration of the poly- 
mer chain. It also depends on the density fluctuation 
length scale fp, which is expressed as a function of Rg 
and Pch as = Rg/[\/2 + 2-KpchR\), and on the pa- 
rameter /lo = h{k = 0) = —(1 — 2(^p/ R^g) / Pch- In fact, 
the dimensionless combination Dj3C,soit is determined by 
only three parameters: p, TV and Rg. In conclusion, once 
thermodynamic parameters, Rg and d are defined, there 
are no adjustable parameters in our method. 

The predicted diffusion coefficients from our rescaled 
MS MD simulations are in good agreement with the data 
for all the test systems. As an example. Table IV displays 



the diffusion coefficients obtained directly from the MS 
MD, D^^^ , once they are rescaled to include the internal 
degrees of freedom, and after the second rescaling of the 
friction, Dcm, as well as the values of the diffusion coeffi- 
cient from the UA MD simulations, D^^, against which 
we compare our predicted diffusion. The table shows that 
while the initial values of the diffusion are orders of mag- 
nitude larger than the data from UA MDs, the rescaled 
coefficients are very close to the real values. For the en- 
tangled systems we adopt the perturbative approach de- 
scribed in Section |VI A| obtaining predicted values that 
are in quantitative agreement with the UA MD simu- 
lations. Entangled samples are from references [T4l I15j . 
and are mostly in the weakly entangled regime. For these 
samples the UA MD simulations include a small number 
of chains: n = 40 for TV = 78, n = 22 for iV = 142, n = 32 
for iV = 174, N = 224, N = 270 and iV = 320, with n 
the number of chains in the simulation and N the degree 
of polymerization. These numbers show one advantage 
of adopting a coarse-grained description as typically our 
samples include thousands of chains. Simulating a large 
ensemble of molecules is necessary, for example, when the 
goal is to investigate large-scale fluctuations or the rel- 
ative relevance of intra- vs inter-molecular contributions 
to the dynamics. 

Fig. [5] illustrates how our approach can be used to cal- 
culate dynamics also in the short time regime. The fig- 
ure shows the mean-square-displacement of the center-of- 
mass from UA MD in comparison with the one calculated 
from the diffusion coefficient rescaled from the MS MD. 
The agreement is quantitative in the long time regime. 
In the short time regime, the UA MD simulation data 
exhibit a subdiffusive behavior, even if polymers are un- 
entangled. In a series of papers we have shown that the 
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TABLE IV: DifTusion coefficients in A^/ns from MS MD com- 
pared with UA MD simulation 



System 


r[Kj 


r,MS 
U 


7") 


L) 


PE 30 


400 


4.44X 10 


58.9 


82.9 


FE 44 


400 


o.iyxio 


44.5 


46.0 


rbj 4o 


450 


o.oOxiO 


36.7 


50.8 


PE 66 


448 


6 n4y 1 0'' 


29.0 


31.8 


PE 78 


450 


6.73x10^ 


23.6 


26.0 


PE 96 


448 


6.98x10^ 


21.9 


23.3 


PE 142 


450 


8.45x10^ 


6.92 


7.93 


PE 174 


450 


8.51x10^ 


4.53 


5.70 


PE 224 


450 


8.80x10^ 


2.73 


3.28 


PE 270 


450 


9.39x10^ 


2.14 


2.06 


PE 320 


450 


8.73x10^ 


1.03 


1.30 



subdiffusive regime is a consequence of the presence of 
cooperative dynamics involving several polymer chains 
moving in a correlated way inside the dynamically het- 
erogeneous hquid of macromolecules.[5^ A detailed 
discussion of this phenomenon, which is of intermolecular 
origin, has been provided before and will not be repeated 
here. Our theory, the Cooperative Dynamics Generalized 
Langevin Equation (CD-GLE), needs as an input the dif- 
fusion coefficient and predicts the subdiffusive behavior 
for times shorter than the longest correlation time as a 
function of the number, n' oc ^/N , of macromolecular 
chains moving in a cooperative way. Fig. [5] shows the 
results from our CD-GLE calculations as dashed lines. 
Here, the input monomer friction coefficient is calculated 
from MS MD using the rescaling procedure. The number 
of chains undergoing cooperative dynamics is n' — 30 for 
N = 96, n' = 25 for N = 66, and n' = 14 for N = 44. 

The subdiffusive behavior shown in UA-MD data is 
not visible in the rescaled MS MD data as the dynamics 
are accelerated. The effective temperature experienced 
by the polymer is much higher than the temperature in 
UA MD simulations as the energy is not dissipated in the 
internal degrees of freedom. 

Finally we discuss the calculations of the monomer and 
soft-colloid friction coefficients for a set of polyethylene 
chains investigated experimentally. [39l I48H5T] The exper- 
imental data do not report the values of Rg at the de- 
sired thermodynamic conditions, T = 509 K and density 
p — 0.0315302 sites/A^, while it is known that the chain 
conformation, and are temperature dependent. To 
calculate the input parameters for our MS MD simula- 
tions we adopt a freely rotating chain model, for which 
the mean-square end-to-end polymer distance is given 
by[67] 



1 + g 2.g 1 - g 



AT 1 



(55) 
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FIG. 5: Plot of mean-square displacement as a function of 
time for unentangled PE melts. The rescaled MS MD simu- 
lation (line) is compared with UA MD simulation (symbols) 
for — 44, 66, 96. Also shown is the outcome of the theory 
for cooperative dynamics (dashed lines). 



For polyethylene melts at this temperature the stiffness 
parameter is g = 0.785. |39j For the samples investi- 
gated in the UA MD simulations, a comparison of the 
theoretical values of Rg using the freely-rotating-chain 
model (FRC) and the values measured directly from the 
UA MD, which are both reported in Table |lj show a rea- 
sonable agreement. The agreement is particularly good 
for the samples that have long chains because for them 
the hypothesis of a Gaussian intramolecular distribution 
is well justified. In this framework, the FRC model pro- 
vides a reliable description of the chain intramolecular 
structure. 

The values of the monomer and soft-colloid friction co- 
efficients for the experimental samples, calculated from 
Eqs. 



(44 1 and (49 1 respectively, are presented in Ta- 
ble |Vj The Table shows the large difference between the 
predicted dimensionless friction coefficients, DjSC^soft and 
NDjiQrm for the same macromolecule coarse-grained at 
two different length scales. From the values displayed in 
Table [V] we calculate the rescaling factor for the friction 
coefficient measured in MS MD simulations, following the 
procedure described in this paper. Because the data have 
different thermodynamic parameters of density and tem- 
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perature, their scaling behavior cannot be inferred from 
their plot, even if an apparent scaling is followed by 
the unentangled samples and the typical reptation N^^ 
scaling by the entangled ones. 



TABLE V: Theoretically calculated dimensionless friction co- 
efficient for monomer (d = 2.1 A) and soft colloid with Rg^^ 
for experimental samples 



System 




D/3Csoft 




PE 36 


101.4350 


0.007846 


0.5153 


PE 72 


219.5710 


0.004927 


0.8946 


PE 106 


331.1465 


0.004543 


1.8407 


PE 130 


409.9056 


0.003497 


1.5354 


PE 143 


452.5669 


0.003318 


1.6822 


PE 192 


613.3669 


0.002828 


2.2412 


PE 242 


777.4485 


0.002500 


2.8186 



T = 509K, p = 0.0315302 [sites/A^^ 



TABLE VL Predicted diffusion coefficients in A^/ns from MS 
MD and experimental data from Refs. |39l I48H51| 



System PE36 PE72 PE106 PE130 PE143 PE192 PE242 
Dcm 111 50 17 14 11 5.6 3.2 
jjexp ^20 41 14 12 8.6 6.5 4.5 



D. Theoretical predictions of diffusion coefficients 
for polyethylene samples 

In this section we report theoretical predictions from 
rescaled mesoscale simulations of the diffusion coefhcients 
for a series of PE samples for which data of chain dynam- 
ics, either from simulations or from experiments, are not 
available in the literature. The degree of polymeriza- 
tion of each sample is not larger than the ones already 
investigated. However, because there are no data to fit 
any parameter, these calculations illustrate the predictive 
power of our method. Diffusion coefficients calculated by 
combining the mesoscale simulations with the rescaling 
procedure presented in this manuscript, are displayed in 
Fig. |6] as a function of the degree of polymerization. 

The set of MS MD simulations is performed for N = 
40,60,80,100,200,300,400,500 and 1,000, at constant 
monomer density pm = 0.0329497 [sites/A^] and tem- 
perature T = 450 K for all samples. Values of Rg are 
calculated using a freely rotating chain model. The hard 
sphere diameter is fixed to the value reported in the pre- 
vious sections for PE, d = 2.1 A, and the pair distribution 
function at contact is g{0) ~ 1/2 as described early on in 
this paper. While the simulations of the small samples 
can be performed on a single CPU machine, for systems 



with a higher degree of polymerization is convenient to 
adopt parallel computing. For those systems, simulations 
were run using the LAMMPS code[7S], with our poten- 
tial as an input, remotely on a 64 CPU machine avail- 
able through the TeraGrid[71]. For the PEIOOO sample, 
which included 85,184 molecules, results were obtained 
after one week of calculations. By comparison with our 
single CPU calculations, running the simulation in par- 
allel reduces the computer time by a factor of 10^. The 
number of particles in the simulation is determined by 
the length of the box size, which for PEIOOO sample is 
equal to 24 Rg, i.e. larger than twice the range of the 
potential, to eliminate molecular self-interaction through 
the periodic boundary conditions. 

While we assume that small changes of density and 
temperature do not affect the hard-core diameter, d, even 
a small difference in p„j can noticeably change the prefac- 
tor in Eq.(27). The monomer friction coefficient is calcu- 



lated using Eq. ( 44 ) for unentangled systems and Eq. ( 54 ) 
for entangled ones. The full lines in Fig. |6] show the 
equations used in the calculation, while the dashed lines 
represents the prediction of the equation for the entan- 
gled system in the unentangled region, and the prediction 
of the unentangled equation in the entangled region. In 
the crossover regime, N w 100 for PE, both expressions 
lead to very similar results, providing a smooth crossover 
between the two equations. 

The predicted values of the diffusion coefficient ap- 
pear to be consistent with the known experimental be- 
havior. The diffusion coefficients of unentangled chains 
{N < 130) follow the scaling behavior of the Rouse ap- 
proach, while the entangled chains show a scaling with 
degree of polymerization of —2.5. Although the latter 
scaling exponent disagrees with the "reptation model", 
it is known that experimental samples of weakly entan- 
gled chains also show a scaling exponent of —2.5. For 
those polymer chains, which are just across the transition 
from unentangled to entangled dynamics, constraint re- 
lease and "tube" fluctuations are relevant. The observed 
scaling behavior is also consistent with the scaling of the 
viscosity observed experimentally. [77J The advantage of 
our method with respect to UA-MD simulations is that 
even in the case of long entangled chains it is possible 
to include a large number of molecules, improving the 
statistics of calculated correlation functions. Overall this 
plot shows that it is possible to provide reasonable pre- 
dictions of large scale dynamical properties by properly 
rescaling mesoscale simulations. 



VII. DISCUSSION AND CONCLUSIONS 

The need for developing a fundamental approach to 
rescale dynamical data obtained from MS simulations of 
coarse-grained systems has been a long-standing problem 
from the time that coarse-graining approaches started be- 
ing developed. Because MS simulations are less compu- 
tational demanding than atomistic simulations, it is pos- 
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FIG. 6: Plot of diffusion coefficients (symbols) as a function 
of degree of polymerization, N. Also shown are the scaling 
exponents for our unentangled, A'^"^ (dotted line), and entan- 
gled systems, N~'^'^ (dotted line). 



sible to investigate larger systems for longer times than 
in all-atom simulations, allowing one to extend the max- 
imum time and length scales accessible through simula- 
tions and to improve the statistics of measured averaged 
quantities. Considering that the number of particles in 
a simulation should be large enough to ensure proximity 
to the thermodynamic limit, MS simulations of coarse- 
grained systems could become an indispensable tool to 
investigate the structure and dynamics of macromolecu- 
lar liquids. 

One advantage of a MS simulation of a coarse-grained 
system is that the simulation speeds up because of the 
averaging of the internal degrees of freedom, leading 
to a softer potential and allowing the study of longer 
timescales than in a fully atomistic simulation. This 
implies, however, that the dynamical properties result- 
ing from the MS MD are faster than their real counter- 
part, for example the ones from UA MD and need to be 
rescaled. 

It is the common procedure to rescale the measured dy- 
namics numerically by bringing a time correlation func- 
tion to agree with the one measured in atomic level sim- 
ulations, however we adopt a different strategy. We have 
proposed a first principle approach to derive an analyti- 
cal form of the rescaling procedure to be applied to the 
dynamics measured directly from MS MD of a, coarse- 
grained polymer liquid. Our approach allows for the re- 
liable prediction of the long-time diffusion of a polymer 
melt as it would bo measured in an atomistic or UA-MD 
simulation. The rescaling procedure has been tested so 
far against simulations and experiments of polyethylene 
liquids both unentangled and entangled. Calculated dif- 
fusion coefficients for samples for which we do not have 
data either experimental or simulated, show consistent 
behavior. 

We start by running MS MD simulations of coarse- 



grained polyethylene melts where each polymer is repre- 
sented as a point particle. The analytical intermolecu- 
lar potential, input to the MS MD, is derived from the 
Ornstein-Zernike equation with the hypernetted closure 
approximation. The correction term to the measured dy- 
namics of the MS MD simulations, is calculated from the 
solution of the Generalized Langevin Equations written 
for the coarse-grained and for the monomer-level repre- 
sentations of the macromolecular liquid. Those equations 
are formally derived from the Liouville equation by as- 
suming two different lengthscales characterizing the rel- 
evant slow dynamics, i.e. monomer and center-of-mass, 
onto which the Liouville equation is projected. 

While the Mori-Zwanzig projection operator technique 
suggests a reliable criteria to select the proper length 
scale of coarse-graining for dynamical properties, the 
GLEs thus generated allow one to derive analytical 
forms of the rescaling contributions associated with the 
coarse-grained dynamical equations. The rescaling pro- 
cedure includes two contributions, given by the changes 
in entropy and in the friction coefficient during coarse- 
graining. The entropic contribution emerges from the 
averaging of the internal degrees of freedom, while the 
friction is due to the change in shape, and as a conse- 
quence the change of the molecular surface exposed to the 
surrounding molecules. Both corrections depend on the 
thermodynamic conditions of the system simulated, and 
on the molecular structure through the radius-of-gyration 
of the iiia.cromokxnile. Thermodynamic and molecular 
quantities enter both directly through the rescaling equa- 
tions and indirectly through the effective potential in the 
mesoscale simulations. In this way the; dynamics pre- 
dicted from the rescaling of each mesoscale simulation is 
specific of the system under study. 

A feature of the coarse-graining models we study is the 

mapping of the polymeric liquid onto simple representa- 
tions, which are isotropic. At the molecular level the 
polymer is described as a soft isotropic sphere. At the 
monomer level, the bead-and spring description affords 
equivalent beads in the chain, which is a reliable approxi- 
mation due to the high number of statistically equivalent 
structural configuration of the molecule. Chain end ef- 
fects enter in the model through the finite size of the 
polymer in the matrix representation of the equations. 
Moreover, because the monomers in a homopolymer are 
structurally identical, with the exception of the two end 
monomers, the intcrmolecular monomer-monomer hard- 
core interaction potential is assumed to be identical for 
any pair of monomers, and each monomer is supposed to 
have identical friction coefficient. 

Although the theoretical picture is straightforward, 
our approach has the advantage of being described in 
closed-form expressions, even if approximated, which al- 
lows for an analytical solution of the rescaling formal- 
ism. This has the potential of being useful in improving 
our understanding of the nature of coarse-graining pro- 
cedures. 
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APPENDIX I: INTERNAL ENERGY 

CALCULATION FOR A 
FREELY-ROTATING-CHAIN MODEL 

The effective mean-force potential for one homopoly- 
mer composed of N monomers can be expressed through 
the structural matrix A as 



2/2 



^,3 



—Y Ar . 

2/2 



(56) 



with the matrix A being real and symmetric, and diag- 
onalized by the orthonormal matrix of the eigenvectors 
Q^i = Q^, so that 



r^Ar = e^Q-i AQC = C^A^ , 



(57) 



where A is the matrix of the eigenvalues, and ^ is the 
matrix of the normal modes defined by r = Q^. 

In this model, the equilibrium distribution function is 



Az 



Ne ^ 



(58) 



where for convenience of notation we introduced the ma- 
trix A' — 3 A/ P. Here Nx is the normalization factor, 
defined by enforcing / dx'^x = 1, as 



/ 3 



V27r/2 



N/2 



[det{A) 



1/2 



(59) 



with Nx ^ Ny ^ ^ N^''^. The statistically aver- 
aged internal energy for one molecule consisting of N 
monomers simplifies to 



kBT/ J J 2 



(60) 



In one dimension, 

E \ _ 



dx X Ax e 



Ax 



= Nx 



2/2 



NxP yr /27r/2 



3A,; 



N 



(61) 



which gives, as the final result for the internal energy of 
one molecule consisting of N monomers, 



E 
kef 



3N 



(62) 



APPENDIX II: THE DYNAMIC MEMORY 
FUNCITION 



We briefly report here the derivation of Eq. ( 29 1 start 



ing from Eq.(13). The product of the direct and pro- 



jected forces is expressed as a function of the density 
field variables as 

(F(0) • FQ{t)) - F(r) • F(r') (p„(r; 0)p^(r; t)) . (63) 

Because the fluid is uniform and isotropic, the den- 
sity fields can be replaced by their fluctuation variables, 
Apct(r,i) = pa{r,t) — (^^(r)), where the ensemble- 
averaged density field is approximated by (pQ.(r)) « 
pg{r). The correlation of the random forces is then ex- 
pressed as 



(F(0) • F'3(<)) ^ f ■ f' p'^g(r)g{r')F{r)F{r') 
(Ap„(r)Ap^(r',0) 

{po.{r)){p,{r')) ' 



(64) 



where we adopt a kind of " dynamical" Kirkwood super- 
position approximation in a weighted average form 

(Ap„(r)Ap^(r',i)) ~ pj dB. g{r)g{r')S{R,t) 

xS-dr-r' + RIi) . (65) 



Eq. ( 65 ) describes the multipoint correlation between the 



density fluctuations at a distance r from segment a at 
time zero, and the density fluctuations a distance r' from 
segment 7 at time t. Because a and 7 can be on the same 
or on different polymer chains, no assumptions are made 
a priori about the relative importance of intra and inter- 
molecular correlations. In this way, the chain connectiv- 
ity does not play a dominant role in our description from 
the very beginning. We then calculate both intra and 
intermolecular contributions and show that both need to 
be included in the calculation of the memory function 
as intramolecular contributions are comparable in size 
to the intermolecular ones. Substitution of Eq.(65l into 
Eq.llell) leads to Eq.([29|. 
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